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ACHILL SCHURMANN 

Abstract. We introduce a parameter space for periodic point sets, given as 
unions of m translates of point lattices. In it we investigate the behavior of the 
sphere packing density function and derive sufficient conditions for local op- 
timality. Using these criteria we prove that perfect, strongly eutactic lattices 
cannot be locally improved to yield a periodic sphere packing with greater den- 
sity. This applies in particular to the densest known lattice sphere packings in 
dimension d < 8 and d = 24. 
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1. Introduction 

The classical and widely studied sphere packing problem asks for a non- 
overlapping arrangement of equally sized spheres in a Euclidean space, such that 
the fraction of space covered by spheres is maximized. The problem arose from the 
arithmetical study of positive definite quadratic forms. By the works Thue [T hulOl 
and Hales |Hal05 1 the optimal arrangements of spheres are known in dimension 2 
and 3. We refer to IIGL87L llCS99l . IIMar03ll and IISch09ll for details and further 
reading. 

For reasons related to the historical roots of the sphere packing problem, special 
attention has been on (point) lattices as the discrete set of sphere centers. In di- 
mension 2 the hexagonal lattice and in dimension 3 the face-centered-cubic lattice 
yield optimal sphere packings. For the restriction of the sphere packing problem 
to lattices, the optimal configurations are known up to dimension 8 and in dimen- 
sion 24 (see Table[T]). Here, solutions are given by fascinating objects, the so-called 
root lattices and the Leech lattice. We refer to liCS99ll . IIMar03ll and llNSl for further 
information on these exceptional objects. 

A major open problem in the theory of sphere packing is to find a dimension 
in which there is a non-lattice packing that is denser than any lattice packing. In 
dimension 10 there exists a non-lattice sphere packing, that is conjectured to have 
a higher density than any lattice sphere packing (see |LS70|). As shown in Table[TJ 
below dimension 24 similar sphere packings have been found in dimensions 11, 
13, 18, 20 and 22. All of them are periodic, that is, a finite union of translates of a 
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lattice sphere packing. By a well-known conjecture, attributed by Gruber HGruO?! 
to Zassenhaus, optimal sphere packing density can always be attained by periodic 
sphere packings. It is known that their density comes arbitrarily close to the optimal 
value (see for example [CE03, Appendix A]). 

A natural idea to obtain a better non-lattice sphere packing, is to "locally mod- 
ify" one of the optimal known lattice sphere packings in dimensions d = 4, . . . , 8. 
In this paper we show that such modifications are not possible within the set of all 
periodic sphere packings (see Corollary [TTI). We more generally show in Theorem 
[TO]that such modifications are not possible for perfect, strongly eutactic lattices. 

One may wonder why the restriction to periodic structures is necessary. One 
could also consider more general discrete sets. However, within the set of all dis- 
crete sets, we are not aware of any notion of a "local modification" that on the one 
hand could potentially lead to an improved sphere packing density, but on the other 
hand would allow us to generalize the result of this paper. For instance, a natural 
approach to define the e-neighborhood of a discrete set is as the collection of sets 
that can be obtained by changing the position of elements by at most an e distance. 
However, such a local modification would not even change the sphere packing 
density. It is equal to a constant multiple of the average number of points per unit 
volume, which could not be changed in such an e-neighborhood. In contrast to 
that, the local changes of periodic sets considered in this paper allow arbitrarily 
large displacements of points, if they are far enough from the origin. 

The paper is organized as follows. In Section |2] we recall some necessary back- 
ground on lattices and positive definite quadratic forms. In Section[3]we introduce 
the so-called Ryshkov polyhedron, and based on it we give a geometrical interpre- 
tation of Voronoi's characterization of locally optimal lattice sphere packings. This 
viewpoint allows a natural generalization to study local optimal periodic sphere 
packings. For their study we introduce a parameter space in Section |4] We give 
characterizations of local optimal periodic sphere packings with up to m lattices 
translates in Section [5] Based on these general characterizations we obtain one of 
the main results of this paper in Section |6l We show that perfect, strongly eutactic 
lattices cannot locally be modified to yield a better periodic sphere packing - they 
are periodic extreme (see Definition [H). 



Lattices and Periodic Sets. A (full rank) lattice L in is a discrete subgroup 
L = Zai + . . . + Ztad generated by d linear independent (column) vectors G 
M^. We say that these vectors form a basis of L and associate it with the matrix 
A = (ai,...,arf) G GLrf(M). We write L = AU^. It is well-known that L 
is generated in this way precisely by the matrices A\J with \J G GL(i(Z). We 
refer to IGL87 1 for details and more background on lattices. Given a lattice L and 
translational vectors ti, for say i = 1, . . . ,m, the discrete set 



2. Background on lattices and quadratic forms 
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d 


point set 




author(s) 


2 


A2 


0.2886... 


Lagrange, 1773, [Lag73| 


3 


As = Ds,* 


0.1767... 


GauB, 1840, [Gau40l 


4 


D4 


0.125 


Korkine & Zolotareff, 1877, fKZ77l 


5 


D5,* 


0.0883... 


Korkine & Zolotai-eff, 1877. I.KZ771 


6 


Ee,* 


0.0721 . . . 


Blichfeldt, 1935, IIBU35II 


7 


E7, * 


0.0625 


Blichfeldt, 1935, [Bli35|| 


8 


Es 


0.0625 


Blichfeldt. 1935. I.Bli35l 


9 


Ag,* 


0.0441 . . . 




10 


PlOc 


0.0390... 


Leech & Sloane. 1970. ||LS70| 


11 


Plla 


0.0351 . . . 


Leech & Sloane, 1970, IILS70I1 


12 


K12 


0.0370... 




13 


PlSa 


0.0351 . . . 


Leech & Sloane, 1970, iLS70i 


14 


Al4,* 


0.0360... 




15 


Al5,* 


0.0441 . . . 




16 


A16,* 


0.0625 




17 


Al7,* 


0.0625 




18 


V18 


0.0750... 


Bierbrauer & Edel, 1998, HBEOOI 


19 


Al9,* 


0.0883... 




20 


V20 


0.1315... 


Vardv. 1995. iVar95l 


21 


A21, * 


0.1767... 




22 


V22 


0.3325 . . . 


Conwav & Sloane. 1996. ICS96l 


23 


A23 


0.5 




24 


A24 


1 


Cohn & Kumar. 2004. IICK09I 



Table 1. Point sets defining best known sphere packings up to dimension 24. In 
dimensions d < 8 and (i = 24 the corresponding authors solved the lattice sphere 
packing problem. The other mentioned authors found the listed, densest known 
periodic sphere packings. The asterisk * indicates that an equally dense, periodic 
non-lattice sphere packing is known. 



is called a periodic (point) set. 

The sphere packing radius A (A) of a discrete set A (not necessarily periodic) 
in the Euclidean space M.'^ (with norm || • ||) is defined as the infimum of half the 
distances between distinct points: 

A(A) = \ inf \\x-y\\. 

The sphere packing radius is the largest possible radius A such that solid spheres 
of radius A and with centers in A do not overlap. Denoting the solid unit sphere by 
B'^, the sphere packing defined by A is the union of non-overlapping spheres 
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Its density (5(A) is, loosely speaking, defined as the fraction of space covered by 
spheres. We can make this definition more precise by considering a cube 

C = {a; G R'' : \xi\ < 1/2} 

and setting 

^(A) = A(A)'^voli?^-liminf ^"^^(t!';^^^ - 

\-^oo vol AC 

If the limit inferior above is a true limit, the cube in the definition can be replaced 
by any other compact set C that is the closure of its interior, without the value of 
5 changing. We say that a corresponding set A is uniformly dense in that case. 
It can be shown that the supremum of 5{h) over all discrete sets is attained by a 
uniformly dense set A. We refer to ||Gro63ll and IICE031 Appendix A] for further 
reading. 

For general discrete sets, it may be difficult to compute the density, respectively 
the limit inferior in the definition. For a lattice the limit inferior can simply be 
replaced by 1/ det L, where det L = \ det A\ is the determinant of the lattice L = 
AU^. Note that the determinant of L is independent of the particular choice of the 
basis A. For periodic sets A as in ([T]) we get the estimate 

,,,, mXikYYoXB'^ 

(5(A) < / , , 

^ ' - det L 

with equaUty if and only if the lattice translates ti + L we. pairwise disjoint. 

Positive definite quadratic forms. Among similarity classes of lattices, hence in 
the space Orf(M)\GLrf(M)/GLrf(Z), there exist only finitely many local maxima 
of (5 up to scaling. In order to characterize and to work with them, i.e., enumerate 
them, it is convenient to use the language of real positive definite quadratic forms 
(PQFs for short). These are simply identified with the set of real symmetric, 
positive definite matrices. Given a matrix Q E 5^q, we set Q[x\ = x^Qx for 
X G W^, defining a corresponding PQF Note that every matrix Q G S^q can be 
decomposed into Q = A* A with A G GLrf(M) and therefore Si,^ can be identified 
with the space Od(M)\GLd(M) of lattice bases up to orthogonal transformations. 
Two PQFs (respectively matrices) Q and Q' are called arithmetically equivalent 
(or integrally equivalent) if there exists a matrix U G GLrf(Z) with Q' = U^QU. 
Thus arithmetical equivalence classes of PQFs are in one-to-one correspondence 
with similarity classes of lattices. 

The arithmetical minimum X{Q) of a PQF Q is defined by 

X{Q) = mill Q[x]. 
xei.''\{o} 

If L = AT,'^ with A G GLrf(M) satisfying Q = A^A is a corresponding lattice, there 
is an immediate relation to the packing radius of L: We have X{Q) = (2A(L))^ 
and therefore 
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where 

^^^^ - (detQ)V^ 

is the so-called Hermite invariant of Q. Note that %{•) is invariant with respect 
to scaling. A classical problem in the arithmetic theory of quadratic forms is the 
determination of the Hermite constant 

T-Ld = sup T-LiQ). 

By the relation described above, it corresponds to determining the supremum of 
possible lattice sphere packing densities. Local maxima of the Hermite invariant 
on S'^Q and corresponding lattices are called extreme. 

3. VORONOI'S CHARACTERIZATION OF EXTREME FORMS 

The Ryshkov polyhedron. Since the Hermite invariant is invariant with respect 
to scaling, a natural approach to maximizing it is to consider all forms with a fixed 
arithmetical minimum, say 1, and minimize the determinant among them. We may 
even relax the condition on the arithmetical minimum and only require that it is at 
least 1. In other words, we have 

■Hd = l/inf(detQ)^/^ 

where 

(2) 7^ = {q G : A(Q) > l} . 

We refer to TZ as Ryshkov polyhedron, as it was Ryshkov | Rys70 1 who noticed that 
this view on Hermite's constant allows a simplified description of Voronoi's theory, 
to be sketched below. 

We denote by S'^ the space of real symmetric matrices, respectively of real qua- 
atic forms in d variables. It is a E 
the usual inner product defined by 



dratic forms in d variables. It is a Euclidean vector space of dimension C^t^) with 



d 

{Q, Q') = I'^'lii = trace(Q • Q'). 
Because of the fundamental identity 

Q[x] = (Q,£C£C*), 

quadratic forms Q € S"^ attaining a fixed value on a given a:; G M'^ \ {0} lie all in 
a hyperplane {affine subspace of co-dimension 1). Thus Ryshkov polyhedra TZ are 
intersections of infinitely many halfspaces: 

(3) n = {QeSiQ: {Q, XX*) > 1 for all x € Z'^ \ {0}}. 

It can be shown that TZ is "locally like a polyhedron", meaning that any intersec- 
tion with a polytope (convex hull of finitely many vertices) is itself a polytope. For 
a proof we refer to HSchOQl Theorem 3. 1]. As a consequence TZ has vertices, edges. 
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facets and in general k-dimensional faces (k-faces). For details on terminology and 
basic properties of polytopes we refer to IIZie97ll . 

Perfect forms. The vertices Q of the Ryshkov polyhedron are called perfect forms. 
Such forms are characterized by the fact that they are determined uniquely by their 
arithmetical minimum (here 1) and its representatives 

MinQ = {x€Z'^: Q[x] = A(Q)}. 

A corresponding lattice is called perfect too. The following proposition due to 
Minkowski implies that the Hermite constant can only be attained among perfect 
forms, i.e., the maximal lattice sphere packing density can only be attained by 
perfect lattices. 

Proposition 1 (Minkowski IIMin05ll ). (detQ)^/'^ is a strictly concave function 
on SyQ. 

For a proof see for example IIGL871 § 39.2]. Note, that in contrast to (det Q)^^^, 
the function det Q is not a concave function on 5>q (see IINel74ll '). However 
Minkowski's theorem implies that the set 

(4) {Q GS^o-detQ>D} 

is strictly convex for D > 0. 

Another property of perfect forms which we use later is the following. 

Proposition 2. If Q ^ S'^ is perfect, then Min Q spans R'^. 

The existence of d linear independent vectors in Min Q for a perfect form Q 
follows from the observation that the rank-1 forms xx^ with x G MinQ have to 
span S'^, since they uniquely determine Q through the linear equations {Q, xx*) = 
\{Q). If however Min Q does not span then these rank-1 forms can maximally 
span a (2) -dimensional subspace of S'^. 

Finiteness up to equivalence. The arithmetical equivalence operation Q 1— )• 
U^QU of GL(f(Z) on S^q leaves X{Q), MinQ and also TZ invariant. In fact, 
GLrf(Z) acts on the sets of faces of a given dimension, thus in particular on the sets 
of vertices, edges and facets of TZ. The following theorem shows that the Ryshkov 
polyhedron TZ contains only finitely many arithmetically inequivalent vertices. By 
Proposition[T]this implies in particular that Tid is actually attained, namely by some 
perfect forms. For a proof we refer to IISch091 Theorem 3.4]. 

Theorem 3 (Voronoi IIVorOTII ). Up to arithmetical equivalence and scaling there 
exist only finitely many perfect forms in a given dimension d>l. 

Thus the classification of perfect forms in a given dimension, respectively the 
enumeration of vertices of the Ryshkov polyhedron up to arithmetical equivalence, 
yields the Hermite constant. Perfect forms have been classified up to dimension 8 
(see i:dSV07II ). 
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Characterization of extreme forms. From dimension 6 onwards not every per- 
fect form is extreme (see MMarOBII '). In order to characterize extreme forms within 
the set of perfect forms the notion of eutaxy is used: A PQF Q is called eutactic 
if its inverse Q^^ is contained in the (relative) interior relint V{Q) of its Voronoi 
domain 

V[Q) = cone{a;a;* : x G MinQ}. 
Here cone M denotes the conic hull 

< aiXi : n G N and Xi G M, Qj > for i = 1, . . . , n > 

of a set M. Note that the Voronoi domain is full-dimensional (of dimension ('^^^)) 
if and only if Q is perfect. Note also that the rank-1 forms xx^ give inequalities 
{Q, xx^) > 1 defining the Ryshkov polyhedron and by this the Voronoi domain of 
Q is equal to the normal cone 

(5) {N eS'^: {N, Q/XiQ)) < {N, Q') for all Q' G U] 

of IZ at its boundary point Q/\{Q). 

Algebraically the eutaxy condition Q^^ G relint V(Q) is equivalent to the exis- 
tence of positive ctx with 

(6) Q'^ = ^ a^xx^. 

x£MmQ 

Thus computationally eutaxy of Q can be tested by solving the linear program 

(7) max amin such that ax > amin and ^ holds. 

The form Q is eutactic if and only if the maximum is greater 0. 

Voronoi IIVor07ll showed that perfection together with eutaxy implies extremal- 
ity and vice versa: 

Theorem 4 (Voronoi IIVor07ll ). A PQF Q G is extreme if and only if Q is 
perfect and eutactic. 

We here give a proof providing a geometrical viewpoint that turns out to be quite 
useful for the intended generalization discussed in the following sections. 

Proof. The function det Q is a positive real valued polynomial on S'^, depending 
on the ('^^^) different coefficients qij of Q. Using the expansion theorem we obtain 

d 

det Q = ^ qf^Qij 

i=l 

for any fixed column index j G {1, . . . , d}. Here, = (—1)*+-' det Qij (with 
Qij the minor matrix of Q, obtained by removing row i and column j) denote the 
coefficients of the adjoint form = (det Q)Q^^ G 5>q of Q. Thus 

(8) grad det Q = (det Q)Q~^ 
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and the tangent hyperplane T in Q of the smooth determinant-det Q-surface 
S = {Q' & S^o ■■ det Q' = det Q} 

is given by 

T = {Q' G S"" : {Q-\Q') = {Q-\Q)}. 
Or in other words, is a normal vector of the tangent plane T of 5 at Q. By 
Proposition [T] and the observation that ^ is convex, we know that S is contained 
in the halfspace 

(9) {Q' gS'' ■.{Q-\Q' -Q)>0}, 

with Q being the unique intersection point of S and T. 

As a consequence, a perfect form Q attains a local minimum of det Q (hence 
is extreme) if and only if the halfspace ^ contains the Ryshkov polyhedron TZ, 
and its boundary meets TZ only in Q. This is easily seen to be equivalent to the 
condition that the normal cone (Voronoi domain) V{Q)ofTZatQ contains Q^^ in 
its interior. □ 

Note that eutaxy alone does not suffice for extremality. However, there exist 
only finitely many eutactic forms in every dimension and they can (in principle) 
be enumerated too (see UMarOSl Section 9.5]). Nevertheless, this seems computa- 
tionally more difficult than the enumeration of perfect forms (see IISto75L ||BM96 I. 
HBatOlll . IIEGS02I ). By the geometry of S and T a eutactic form attains always a 
unique minimum of 6 (maximum of det) on its face of the Ryshkov polyhedron. 
However, not all faces of the Ryshkov polyhedron contain a eutactic form. 

4. Parameter spaces for periodic sets 

We want to study the more general situation of periodic sphere packings. Recall 
from (dl) that a periodic set with m lattice translates (an m-periodic set) in M'^ is of 
the form 

m 

(10) A' = U {t'i + L) , 

i=l 

with a lattice L dW^ and translation vectors t'- G W^, i = 1, . . . ,m. 

We want to work with a parameter space for m-periodic sets similar to S^q for 
lattices. For this, we consider A' as a linear image A' = AAt of a standard periodic 
set 

m 

(11) Kt = \j(t,+z'^y 

i=l 

Here, A G GLrf(M) satisfies in particular L = AU^. Since we are only interested 
in properties of periodic sets up to isometrics, we encode A' by Q = A* A S S^q, 
together with the m translation vectors ti, . . . , tm- Since every property of periodic 
sets we deal with here is invariant up to translations, we may assume without loss 
of generality that tm = 0. Thus we consider the parameter space 

(12) 5^'o"' = 5^o xlK'^''^"'"^^ 
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for m-periodic sets (up to isometries). We hereby in particular generalize the space 
= SyQ in a natural way. We call the elements of 5^'^ periodic forms and 
denote them usually hy X = {Q, t), where Q € S^q and 

t = (ti,...,t^_i) GM'^^^^"-!) 

is a real valued matrix containing m—1 columns with vectors ti G M.'^. One should 
keep in mind: although we omit tm = 0, we implicitly keep it as a translation 
vector. Note that a periodic set A' as in (flOl ) has many representations by periodic 
forms. In particular, m may vary and we have different choices for A. A similar 
approach for periodic sets in dimension 3 has been considered in IIPZ98II . 
The parameter space 5^'™ is contained in the space 

(13) gd,m ^gd ^^dx(ni-l)_ 

It can be turned into a Euclidean space with inner product (•,•), defined for X = 
{Q,t) and a:' = (Q',t') by 

m—l 

{X,X') = {Q,Q') + Y,t% 
1=1 

Note, for the sake of simplicity we use the same symbol for the inner products on 
all spaces 5'^'™. 

We extend the definition of the arithmetical minimum A, by defining the gener- 
alized arithmetical minimum 

\{X) = min{Q[ti - tj - v] : I < j < m and v ^ Z'^, with v ^Oiii= j} 

for the periodic form X = {Q,t) G 5^'^. Note that we have X{X) = in the 
case of intersecting lattice translates (tj + Z'^) n (tj + Z'^) ^ with i ^ j. The 
set of representations of the generalized arithmetical minimum Min X is the set 
of all w = ti — tj — V attaining X{X). Computationally, MinX and XiX) can 
be obtained by solving a sequence of closest vector problems (CVPs), one for each 
pair i, j with i ^ j. In addition one shortest vector problem (SVP) has to be solved, 
taking care of the cases where i = j. Implementations of algorithms solving CVPs 
and SVPs are provided for example in MAGMA [MAG] or GAP [ GAP] . 

In order to define the sphere packing density function 6 : 5^'^ — )• M we set 
det X = det Q for periodic forms X = {Q,t). Then 

In analogy to the lattice case, we call a periodic form X G m-extreme if it 

attains a local maximum of 6 within 5^'™. 

The relation (fT4l ) shows that the supremum of 5 among m-periodic sphere pack- 
ings is up to some power and a constant factor equal to the "Hermite like constant" 



sup 



X{X)/{detXf/'^ = 1/ inf (detX)^/'^, 
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where the set TZm on the right side is the (generalized) Ryshkov set 

(15) = |x G S^;;" : A(X) > 1} . 

The condition \{X) > 1 gives infinitely many linear inequalities 

p^{X)=Q[v] = {X, {vv\0)) > 1 

for V G {0}, as in the case m = 1. For m > 1 we additionally have the 
infinitely many polynomial inequalities 

(16) Pij>(X) = Q[i, -tj-t;] > 1, 

where j G {1, . . . , m} with i 7^ j and v € Z'^. These polynomials are of degree 3 
in the parameters qki, tki of X. Note that they are linear for a fixed t. Observe also 
that pi^m,v and Pm,j,v are special due to our assumption tm = and that there is 
a symmetry Pij^v = Pj,i,-v by which we may restrict our attention to polynomials 
with i < j. For i = j we have the linear function Pij^v = Pv 

5. Local analysis of periodic sphere packings 

Characterizing local optima. Before we generalize perfection and eutaxy to a 
notion of m-perfection and m-eutaxy (in order to obtain a sufficient condition for 
a periodic form to be m-extreme from it) we discuss a rather general setting: As- 
sume we want to minimize a smooth function on a basic closed semialgebraic set, 
that is, on a region which is described by finitely many (non-strict) polynomial in- 
equalities. Let E denote a Euclidean space with inner product (•,•). Further, let 
/ : ii^ — )• M be smooth (infinitely differentiable) and gi, . . . ,gj. be (real valued) 
polynomials on E. Assume we want to determine whether or not we have a local 
minimum of / at Xq on the boundary of 

(17) G = {X (^E : gi{X) > for i = 1, . . . , A;}. 

For simplicity, we further assume (grad/)(Xo) ^ and 5i(Xo) = 0, as well 
as (grad 5j)(Xo) / 0, for i = 1, . . . , /c. Then, in a sufficiently small neighborhood 
of Xq, the function / as well as the polynomials gi can be approximated arbitrarily 
close by corresponding affine functions. For example, / is approximated by the 
beginning of its Taylor series 

/(Xo) + ((grad/)(Xo),X-Xo). 

From this one easily derives the following well-known criterion (see for example 
IIBSS93I Theorem 4.2.2]) for an isolated local minimum of / at Xq, depending on 
the normal cone 

V(Xo) = cone{(grad5i)(Xo) : i = 1, . . . , A;}. 
The function / attains an isolated local minimum on G if 

(18) (grad/)(Xo) GintV(Xo), 
and / does not attain a local minimum if 

(19) (grad/)(Xo) 0V(Xo). 
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The behavior in the case (grad/)(Xo) € bdcone V(Xo) depends on the involved 
functions / and gi and has to be treated depending on the specific problem. 

For the lattice sphere packing problem we have E = S'^ and / = det^'''^. For 
<5o e we set gi{Q) = Q[vi] - \{Qq) with {gTa.dgi){Q) = Viv\ for each pair 
±.Vi in MinQo- By Theorem|4]we have a local minimum of f{Q) = (det QY^'^ 
at Qq on G (as in ([17] )) if and only if Qq is perfect and eutactic, respectively if 
V(Qo) is full-dimensional and (grad/)((5o) £ int V((5o)- Here, (grad/)(Qo) is 
a positive multiple of Qq^- Thus in this special case (due to Proposition [T]) we do 
not have a local minimum of / where (grad f){Qo) G bd cone V{Qo). 

Let us consider the case of m-periodic sets, hence of = S'^'"^ with m > 1. 
We want to know if a periodic form Xq G 5^'^ attains a local minimum of / = 
det^/*^. We may assume X{Xo) > 0. The set MinXo is finite and moreover, for 
X = {Q, t) in a small neighborhood of Xq = {Qq, t"), every ti — tj — v G Min X 
corresponds to a. — — v G MinXo. Thus locally at Xq, the generalized 
Ryshkov set TZm is given by the basic closed semialgebraic set G defined by the 
inequalities pij,y{X) — \{Xq) > 0, one for each pair ±(t? — — v) in MinXo. 
As explained in Section IH we may assume 1 < i < j < m and = if j = m. 
An elementary calculation yields 

(20) (gradpij- „)(X) = {ww\ 0, . . . , 0, 2Qw, 0, . . . , 0, -2Qw, 0, . . . , 0), 

where we set X = {Q, t) and use w to abbreviate tj — tj — v. This is to be under- 
stood as a vector in S'^'"^ = S'^ x R^^C*"-!), with its "^'^-component" being the 
rank-1 form ww'^ and its "translational-component" containing the zero-vector 
in all but the ith and jth column. If j = m, the jth column is omitted and if i = j 
the corresponding column is 0. For (grad/)(X) we obtain a positive multiple of 

{Q-\o). 

A sufficient condition for local m-periodic sphere packing optima. General- 
izing the notion of perfection, we say a periodic form X = {Q, t) G 5^'™ (and 
a corresponding periodic set represented by X) is m-perfect if the generalized 
Voronoi domain 

(21) V{X) = cone{(gradpij-^)(X) : ti - tj - v € MinX for some v G Z'^} 

is full-dimensional, that is, if dim V(X) = dim 5"'''^ = (''+^) + (m - l)d. Gen- 
eralizing the notion of eutaxy, we say that X (and a corresponding periodic set) is 
m-eutactic if 

(Q"\0) G relintV(X). 

So the general discussion at the beginning of this section yields the following suf- 
ficient condition for a periodic form X to be isolated m-extreme, that is, for X 
having the property that any sufficiently small change which preserves \{X), nec- 
essarily lowers 5{X). 

Theorem 5, If a periodic form X G S'^^ is m-perfect and m-eutactic, then X is 
isolated m-extreme. 
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Note that the theorem gives a computational tool to certify isolated m- 
extremeness of a given periodic form X = {Q,t) £ 5^'™: First, we compute 
Min X and use equation (l20l ) to obtain generators of the generalized Voronoi do- 
main V{X). From the generators it can be easily checked if the domain is full- 
dimensional, hence if X is m-perfect. Next, we can computationally test whether 
{Q^^, 0) is in V{X) or not; for example by solving a linear program similar to (|7]l. 

If we find G relintV(X) (or equivalently in intV(X) as V(X) is 

assumed to be full-dimensional), the periodic form X represents an isolated m- 
extreme periodic set. If {Q^^, 0) V{X), the periodic form X does not represent 
an m-extreme periodic set. In this situation, we can even find a "direction" N £ 
gd,m^ for which we can improve the sphere packing density of the periodic form X, 
that is, such that 5{X + eN) > 5{X) for all sufficiently small e > 0. 

Remark 6. Let X £ S^^ with V(X). Then we can improve the 

sphere packing density of X in direction N given by the nearest point to —{Q^^ ,0) 
in the polyhedral cone 

(22) V{X) = {N £ 5'^'™ : {V,N) > Ofor ally £ V{X)}. 

Note that the cone V{X) is dual to the generalized Voronoi domain V{X) and 
(added to X) gives locally a linear approximation of the generalized Ryshkov 
set IZm- 

Fluid diamond packings. For general m we are confronted with a difficulty 
which does not show up in the lattice case m = \: There may be non-isolated 
m-extreme sets, which are not m-perfect. The. fluid diamond packings in dimen- 
sion 9, described by Conway and Sloane in [CS95J, give such an example. 

Example. The root lattice can be defined by 

d 

Drf = {:e G : ^ Xi = mod 2}. 

2 = 1 

The fluid diamond packings are 2-periodic sets 

D9(t) = D9U(D9+t) 

with t G such that the minimal distance among elements is equal to the 
minimum distance of Dg itself. We may choose for example t = ta = 
. . . , ^, a)* with any a G M. For integers a we obtain the densest known pack- 
ing lattice Ag = ^^{ta) in dimension 9, showing that it is part of a family of 
uncountably many equally dense 2-periodic sets. 

The sets Dg(ta) give examples of non-isolated 2-extreme sets, which are 2- 
eutactic, but not 2-perfect. In order to see this, let us consider a representation 
Xa £ for D<j(ta). We choose a basis A of Dg. Then Xa = {Q, A~^ta), with 
Q = A^A, is a representation of Dg{ta)- 

For non-integral a we find Min Xa = Min Q (using MAGMA for example). It 
follows (for example by Lemma |9] below) that Xa is 2-eutactic, but not 2-perfect. 
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MinXa = MinQ U {(xi, . . . ,X8, 0)* G {0, 1}^ : ^ 



Xi = mod 2}. 



i=l 



Thus the vectors in Min Xa \ Min Q span only an 8-dimensional space. Therefore 
Xa is not 2-perfect. Nevertheless, a corresponding calculation shows that X^. is 
2-eutactic, as in the case of non-integral a. 

In order to see that X^. is non-isolated 2-extreme, we can apply Proposition |7] 
below. One easily checks that for integral a (hence for the lattice Ag) we have 
only one degree of freedom for a local change of ta giving an equally dense sphere 
packing. For non-integral a we have nine degrees of freedom for such a modifica- 
tion. 



Non-isolated m-extreme sets as in this example can occur for periodic 
forms X G only if {Q^^,0) G bdV(X) (which is for example always 

the case if X is m-eutactic, but not m-perfect). In this case it is in general not 
clear what an infinitesimal change of X in a direction N G S"^'"^ leads to (already 
assuming it is orthogonal to {Q~^, 0) as well as in the boundary of the set V{X) 
in (|22)) ). If T{X) denotes the unique face of V{X) containing {Q^^, 0) in its rel- 
ative interior, then this "set of uncertainty" is equal to the face of V{X) dual to 
J^{X), that is, equal to 



Or in other words, the set U{X) is the intersection of V{X) with the hyperplane 
orthogonal to {Q^^ , 0). Note that it is possible to determine F{X) (and hence a 
description ofU{X) by linear inequalities) computationally, using linear program- 
ming techniques. 

Purely translational changes. Below we give an additional sufficient condition 
for m-extremeness. For this we consider the case when all directions in U{X) 
are ''purely translational changes" N = (0, t^) G 5^^'*". A vivid interpretation 
of a purely translational change can be given by thinking of the corresponding 
modification of a periodic sphere packing. The spheres of each lattice translate are 
jointly moved. If in such a local change all contacts among spheres are lost, we 
can increase their radius and obtain a new sphere packing with larger density. If 
some contacts among spheres are preserved however, the sphere packing density 
remains the same. The latter case is captured in the following proposition, which 
gives an easily testable criterion for m-extremeness. We apply this proposition in 
Section [6l where we consider potential local improvements of best known packing 
lattices to periodic non-lattice sets. 

Proposition 7. For a periodic form X = {Q,t) G 5^'™ with {Q^^,0) G 
hdV{X), letU(X) be contained in 

{(0,t^) G 5"'''" : tf = for at least one ti - tj - v e MinX with v G Z'^}. 
Then X is (possibly non-isolated) m-extreme. 



□ 



(23) 



U{X) = {N e V{X) ■ (F, iV) = for all V G F{X)]. 
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Note, if X is m-eutactic (possibly not m-perfect), the set U{X) is the orthog- 
onal complement V{X)^ of the hnear hull of V{X). Note also that Proposition |7] 
includes in particular the special case where some € Z'^ are in Min X (and there- 
fore ti = tj = for i = j = m). This situation occurs for the 2-periodic, fluid 
diamond packings in the example above. 

From the sphere packing interpretation of the proposition its assertion is clear. 
Nevertheless, we give a proof below, based on a local analysis in S'^]^ . More 
than actually needed for the proof, we analyze how 5 changes locally at a periodic 
form X G 5!^'™ in a direction N G U{X). As a byproduct, we obtain tools 
allowing a computational analysis of possible local optimality for a given periodic 
form (not necessarily covered by the proposition). These can for example be used 
in a numerical search for good periodic sphere packings. 

Proof of Proposition^ The generalized Voronoi domain V{X) is spanned by gra- 
dients (gradpjj^„)(X) (as given in (l20l)). one for each pair of vectors ztw G 
MinX. The assumption that a direction N = {Q^ , t'^) is in U{X) for a periodic 
form X = {Q, t), implies {Q"^ , Q^) = 0. Moreover, for the unique maximal face 
F{X) of V{X) with 0) G relint the condition that N is orthogonal 

to some (gradpjj-^„)(X) in F{X) translates into 

(24) ((gradpij>)(X), N) = Q^H + 2(tf - tffQw = 0, 

with w = {ti — tj — v). Recall that in the special case i = j (and for m = 1 
anyway) Pij,v is linear and (l24l ) reduces to the condition Q^[w] = 0; if then N 
satisfies this linear condition, Pi j> {X + eA^) is a constant function in e. 

When pij^v{X + eA^) is a cubic polynomial in e we need to use higher order 
information in order to judge its behavior. An elementary calculation yields for the 
Hessian 

(25) (hessp,,,, „)(X)[iV] = 2Q[tf - tf ] + 4(tf - iffQ^'w. 

Now how does 6 change at X in direction A^, assuming it is in the set of un- 
certainty U{X)1 Among the polynomials with satisfying (|24l) . the fastest 
decreasing polynomial in direction A^ determines \{X + eN) for small enough e. 
Thus for the local change of 6 in direction A^, we may restrict our attention to a 
polynomial Pij^v with the smallest value (1251 ) of its Hessian. 

By Proposition [T] we know that det^^'^ decreases strictly at X in a direction A^ G 
UiX) if and only if / 0. 

For a purely translational change with = 0, the function det^/*^ remains 
constant. On the other hand, because of (1251 ) and since Q is positive definite, we 
have {hess Pi j^v)iX)[N] > 0, with equality if and only if tf - 1 = 0. The latter 
implies that pij^y{X + eN) is a constant function of e. Thus for purely translational 
changes A^ = (0,t^) G U{X), the density function 6{X + eA^) is constant for 
small enough e > 0, if tf' = tj^ for some pair with ti — tj — v G Min X (for 
a suitable v G Z'^). This proves the proposition. □ 
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Note that our argumentation in the proof also shows that 6{X + eN) increases 
for small e > 0, for a purely translational change N = (0,t^) G ^{X) with 
tf ^ for all pairs with U - tj - v G MinX (for some v G Z'^). 
This case corresponds to a modification of a periodic sphere packing in which all 
contacts among spheres are lost. 

6. Periodic extreme sets 

A given periodic set has many representations by periodic forms, in spaces 5^'™ 
with varying m. For example, by choosing some sublattice of Z'^, we can add 
additional translational parts. 

It could happen that a periodic set A with a given representation X G 5^'^ is 
m-extreme, whereas a second representation X' G 5^^'™' is not m'-extreme. We 
are not aware of an example though. However, in some cases we are certain that 
the packing density of no representation of A can locally be improved. 

Definition 8. A periodic set is periodic extreme if it is m-extreme for all possible 
representations X G S^]^. 

Theorem [TO] below gives a sufficient condition for a lattice to be periodic ex- 
treme. For its statement we need the notion of strong eutaxy for lattices, respec- 
tively PQFs: A form Q G Si,^ (and a corresponding lattice) is called strongly 
eutactic if 

(26) Q'^ = a XX* 

x&A'mQ 

for some q > 0, i.e., if the coefficients in the eutaxy condition ^ are all equal. 
It is well-known that a PQF Q is strongly eutactic if and only if the vectors in 
Min Q form a so-called spherical 2-design with respect to Q (see [.VenOlJ . l,Mar03l 
Corollary 16.1.3]). 

Lemma 9. Any representation X G 5^'^ of a strongly eutactic lattice ( respectively 
PQF) is m-eutactic. 

Proof. Let Q G 5>q be strongly eutactic, satisfying (l26l ) for some a > 0. Let 
X = {Q-^ ,t^) G 5^'™ be some representation of Q, e.g. with m > 1. Let the 
corresponding eutactic lattice be denoted by A. Then is the Gram matrix of a 
basis A G GLrf(M) of a sublattice L of A. The columns of are the coordinates 
of lattice points of A relative to A. 

For a fixed w G MinX we define an abstract graph, whose vertices are the 
indices in {1, ... , m}. Two vertices i and j are connected by an edge whenever 
there is some v such that w = tf —t^ — v. In other words, the graph reflects 
via an edge that spheres of packing radius A(A) around points of the two 
sublattice ti-anslates A(tf + Z'^) and A(tf + Z'^) touch. For z G Z'^, the sphere 
with center A{t^ + z) touches the sphere with center A{t'j- + z + w). Since 
the periodic form X represents a lattice A, we find a chain of touching spheres 
at centers ^(t^" + z + kw), with k = 0, 1, . . .. Modulo some natural number 
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less or equal to m these centers belong to the same lattice translate of L. As a 
consequence, we find that the graph defined above is a disjoint union of cycles. So 
w induces a partition (/i, . . . , Ik) of {1, ... , m}. 

Let / be an index set of this partition (containing the indices of a fixed cycle of 
the defined graph). Summing over all triples {i,j,v) with and G such 

that w = tf -tf -V £ Min X, we find (using (lIOl)): 

^ (gradpi,,>)(X) = 2|I|(W,0). 

{i,j,-u)SI^xZ'' 
with -u—t-S —t^ ~w 

The factor 2 comes from the symmetry gradp, j .„ = gradpj j Summation over 
all index sets / of the partition yields 

(27) (gradpij,„)(X) =2m(W,0). 

(i,J,i')6{l,---,™}^ XZ"* 
with u=t?f -t-?^ —w 

As a consequence we find by the strong eutaxy condition (l26l) that 

{Q~\Q) = (a/2m) (grad^j, 

mSMin J>C,(i ,j,-u)S{l,...,m}^ XZ'' 
with ■!;=t<^-t-''^--!i) 

with a suitable a > 0. Thus X is m-eutactic. □ 

Not all PQFs (or lattices) which are strongly eutactic have to be perfect. For 
example the lattices 17^ forn > 2 are of this kind. But if a strongly eutactic PQF is 
in addition also perfect, then the following theorem shows that this is sufficient for 
it to be periodic extreme. Note that this applies in particular to so called strongly 
perfect lattices and PQFs. For these lattices the vectors in Min Q form a spherical 
4-design with respect to Q (see IINeb02l or llMar03[ Chapter 16] for further details). 

Theorem 10. Perfect, strongly eutactic lattices (respectively PQFs) are periodic 
extreme. 

Proof. Let Q G be perfect and strongly eutactic. Hence the vectors in Min Q 
span (by Proposition [J]) and satisfy (l26]l for some a > 0. Let X = {Q^,t-^) e 
^S^^'™ be a representation of Q. By Lemma |9l X is m-eutactic. If X is m-perfect 
as well, we know by Theorem [5] that X is also m-extreme. 

So let us assume that X is not m-perfect; hence the generalized Voronoi do- 
main V{X) is not full-dimensional. We want to apply Proposition |7j For this we 
choose 

AT = (Q^, t^) G U{X) = V{X)^ with / 0. 

(Recall the definition of U{X) from ^ and that U{X) = V(X)-L if X is m- 
eutactic.) By this assumption we have in particular 



{N, (gradpi,,>)(X)) = 
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for all triples {i,j,v) with w = tf — — v £ MinX. Using equation (|27] ). 
which we obtained in the proof of Lemma |9l we get {N, {ww\ 0)) = Q^[w] = 
for every fixed w G Min X. 

By Proposition |2] there exist d linearly independent w in MinX, which implies 
= 0. Using dm), we obtain 

(28) = {N, (gradpij- = 2(tf - tfYQw. 

If — tj^ = for some pair (i, j) we can apply Proposition |7J Note that this 

includes in particular the case i = j = m {tf = tf = 0) if G Z*^ n MinX. So 
we may assume that such v do not exist. 

We choose d linearly independent vectors wi, . . . , € MinX (that exist by 
Proposition|2l). By the assumption that non of the Wi is integral and by the assump- 
tion that X represents a lattice, each Wi connects the origin t;^ = to another 
translation vector (with j ^ m) via Wi = — v for some v G Z^. In the 
same way each of the chosen minimal vectors connects the translation vector with 
index i to other translation vectors. We denote by / the subset of {1, ... , m} that 
is connected to the index m (respectively to the origin 0) via a sequence of such 
links through the chosen d minimal vectors. For each index i ^ I we get from the 
minimal vectors d independent linear conditions (l28l) for the differences tf — , 
with suitable j ^ I \ {i}. Overall we obtain d\I\ independent equations for d\I\ 
differences. We deduce that all of them vanish. Moreover, as t = we even find 
tf = for all indices i G /. 

□ 

The root lattices Ad, and E^, as well as the Leech lattice are known to be 
perfect and strongly eu tactic (cf. [Mar03]). These lattices are known to solve the 
lattice sphere packing problem in dimensions d < 8 and d = 2A (see Table [T]l- 
As an immediate consequence of Theorem [TOl we find that they cannot locally be 
improved to a periodic non-lattice set with greater sphere packing density. 

Corollary 11. The lattices A^,, for d > 2, D^, far d > 3, and E^, far d = 6, 7, 8, 
as well as the Leech lattice are periodic extreme. 

We also checked whether or not Theorem [TO] can be applied to other dimen- 
sions d < 24. For these dimensions the so-called laminated lattices Ad and sec- 
tions Kd of the Leech lattice give the densest known lattice sphere packings. The 
lattices Kd are different from (and at the same time give the densest known 
lattice sphere packings) only in dimensions (i = 11, 12, 13. For these d, the lattice 
Kd is strongly eutactic only for d = 12, when Kd is also known as Coxeter-Todd 
lattice. The laminated lattices A^ give the densest known packing lattices in di- 
mensions d = 9, 10 and d = 14, . . . , 24 (for d = 18, . . . , 24 they coincide with 
Kd). Among those values for d, the laminated lattices A^ are strongly eutactic 
if and only if d = 15, 16 or d > 20. Concluding, we cannot exclude that densest 
known lattice sphere packings in dimensions d G {9, 10, 11, 13, 14, 17, 18, 19} can 
locally be improved to better periodic sphere packings. Further analysis is required 
here. 
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7. Floating and strict periodic extreme lattices 

The last step of the proof of Theorem [10] has a vivid interpretation if we think 
of a sphere packing described by the given lattice. Let X = {Q-^,t^) be one 
of its representations and let A denote a sublattice basis with Gram matrix . 
Then the sublattice translates A{tf + Z'^) with i £ I form a "rigid component" of 
the sphere packing. If we do not want to decrease the sphere packing density in a 
local deformation we have to move all of its translates simultaneously. This rigid 
component may actually be larger than the one used in the proof of Theorem[lOl It 
may even consist of the whole packing. A maximal rigid component of translates 
can be described via an abstract graph with vertices in {1, ... , m}: is an edge 
whenever there is some v Z'^ such that — — v £ Min X. Let / be the 
set of indices i (vertices of the graph) connected by a path with m. If |/| = m the 
whole packing forms one rigid component. If / is a strict subset of {1, ... , m} we 
say a corresponding packing or lattice is floating. 

In a floating packing each connected component of the graph above corresponds 
to a union of translates which can jointly locally be moved without changing X{X) 
and 6{X) respectively. Examples are the fluid diamond packings described in the 
example of Section[5] The same applies to their higher-dimensional generalizations 
D+ = Dd U (Drf + (i, . . . , i)) for d > 10 (see IICS99. Section 4.7.3]). For even d 
these 2-periodic sets are actually lattices (hence 1-periodic). In fact Eg = Dg. 

We note that Theorem [10] and Corollary [TT] give statements about local opti- 
mality of lattices, but not about strict local optimality. With the assumptions of 
Theorem [lO] alone strict local optimality cannot be ensured, as shown by floating 
lattices like Dj" for even d > 10. These lattices have the same minimal vectors 
as the corresponding root lattice and therefore they give a series of perfect and 
strongly eutactic lattices that are periodic extreme by Theorem [TO] However, they 
can locally be modified to other 2-periodic sets of the same density. 

We think that a strengthening of Corollary [His possible for certain lattices that 
are non-floating, perfect and strongly eutactic. These include in particular the Eg 
root lattice and the Leech lattice. We think it is possible to show that these lattices 
are strict periodic extreme, meaning they are isolated m-extreme for all possible 
representations X G ^S^^,^. By Lemma [9] and Theorem [5] one has to show that a 
given non-floating, perfect and strongly eutactic lattice is m-perfect for every m. 
Here some further work is required... 
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